Atom Làsers, Coherent States, and Coherence: 
II. Maximally Robust Ensembles of Pure States 
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As discussed in the preceding paper [Wiseman and Vaccaro, |quant-ph/9906125 ] , the stationary 
state of an optical or atom làser far above threshold is a mixture of coherent field states with 
random phase, or, equivalently, a Poissonian mixture of number states. We are interested in which, 
if either, of these descriptions of p ss as a stationary ensemble of pure states, is more natural. In 
the preceding paper we concentrated upon the question of whether descriptions such as these are 
physically realizable (PR). In this paper we investigate another relevant aspect of these ensembles, 
their robustness. A robust ensemble is one for which the pure states that comprise it survive 
relatively unchanged for a long time under the system evolution. We determine numerically the 
most robust ensembles as a function of the parameters in the làser model: the self-energy \ °f the 
bosons in the làser mode, and the excess phase noise v. We find that these most robust ensembles 
are PR ensembles, or similar to PR ensembles, for all vàlues of these parameters. In the ideal 
làser limit (y = \ = 0), the most robust states are coherent states. As the phase noise or phase 
dispersion is increased through v or the self-interaction of the bosons Xi respectively, the most 
robust states become more and more amplitude-squeezed. We find scaling laws for these states, 
and give analytical derivations for them. As the phase diffusion or dispersion becomes so large that 
the làser output is no longer quantum coherent, the most robust states become so squeezed that 
they cease to have a well-defined coherent amplitude. That is, the quantum coherence of the làser 
output is manifest in the most robust PR ensemble being an ensemble of states with a well-defined 
coherent amplitude. This lends support to our approach of regarding robust PR ensembles as the 
most natural description of the state of the làser mode. It also has interesting implications for atom 
làsers in particular, for which phase dispersion due to self-interactions is expected to be large. 

03.65.Yz, 03.75.Fi, 42.50.Lc, 03.65.Ta 
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I. INTRODUCTION 

A làser is a device that produces a coherent beam of 
bosons. The meaning of the word 'coherent' in this con- 
text is discussed at length in a paper by one of us |Q. 
In particular, a coherent output does not mean that the 
output, or the làser mode itself, is in a coherent state. 
Rather, as has long been recognized 0, the stationary 
state matrix for the làser mode is a mixture of number 
states. In the far-above threshold limit, this mixture is 
Poissonian with mean ji: 
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(1.1) 



This state matrix can also be represented as a mixture 
of coherent states: 



2tt 



><l< 



a e a e 



I 



(1.2) 



where |a| 2 = fx. 



On the basis of this second representation, one might 
claim that the làser really is in a coherent state | |a|e l( ^), 
but that one cannot know a priori what the phase (j> is. 
In the preceding paper || we have investigated whether 
this claim is true. If it were true then there should be 
some way of finding out which coherent state the làser 
is in without affecting its dynamics. We found that if 
there is any self-energy in the làser mode (such as a x 
nonlinearity for an optical làser, or s-wave scattering for 
an atom làser), then it is in fact not possible tophysi- 
cally realize the coherent state ensemble in Eq. (1.2). By 
contrast, it is always possible to physically realize the 
number state ensemble in Eq. (1.1). 



For an ideal làser (with no x^-like nonlinearity), the 
unknown coherent state description and the unknown 
number state description are both physically realizable 
(PR). Given that they arc mathcmatically equivalent, 
why is the former description ubiquitous and the latter 
rare? The answer, as was pointed out some time ago 
by Gea-Banacloche M, is differential survival times. An 
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ideal làser prepared in a coherent state will remain close 
to that initial state for a time of order K , where k is 
the bare decay rate of the cavity. By contrast, a làser 
prepared in a number state will be likely to remain in 
that state only for a time of order where \i is the 

mean number as above. 

This result, derived also in Ref. Q, was taken fur- 
ther by Gea-Banacloche in Ref. ||] using the early model 
for a làser with saturation due to Scully and Lamb 
Gea-Banacloche considered pure states with mean pho- 
ton number equal to that of the làser at steady state, and 
calculated their purity at later times. He showed that the 
pure state that had the slowest initial rate of decay of pu- 
rity was, in general, a slightly amplitude-squeezed state 
rather than a coherent state. 

There seems little doubt, then, that it is most useful 
to consider an ideal làser to be in a coherent state (or 
nearly coherent state) of unknown phase. However it is 
an open question whcthcr this is true of a non-ideal làser, 
that is, a làser with additional noise or dispersion of some 
form. Another open question is how this issue relates to 
the quantum coherence of the output of such a non-ideal 
device. 

The particular làser system of interest here is the atom 
làser jï|. An important difference between an atom làser 
and an optical làser is that the interatomic interactions 
cannot be neglected. This gives rise to a x^-like nonlin- 
earity in the làser mode. As noted above, this affects the 
physical realizability of ensembles, and we also expect it 
to affect their robustness. 

A robustness analysis for a Bose-Einstein condensate 
has been done by one of us with Barnett and Burnett . 
This produced similar results to that of Gea-Banacloche 
||, although it was based on the fidelity || which mea- 
sures the overlap of the initial state with the state at a 
later time. However, the authors of Ref. only calcu- 
lated the initial rate of decay of the fidelity, and this is 
unaffected by any Hamiltonian terms. Hence the self en- 
ergy played no role in this analysis. Moreover, the treat- 
ment, like that of Gea-Banacloche 0, considered only a 
single pure state to represent the state of the condensate. 
Thus it does not give, in gene ral, a represe nta tion of the 
steady state on par with Eq. (El) or Eq. (1.2). 

In this paper we give an analysis that treats the dy- 
namics of an atom làser at all times and that incorporates 
an ensemble of pure states. It takes into account Hamil- 
tonian terms and gives a robust representation of the 
steady state. We consider both the problcm of finding 
the most robust ensemble, and the most robust physi- 
cally realizable (PR) ensemble. Since ensembles are re- 
alized by unraveling the master equation finding 
the most robust PR ensemble is equivalent to finding the 
maximally robust unraveling, a concept introduced by us 
in Ref. §. 

A review of maximally robust unravelings, including a 
comparison with other approaches, is given in Sec. II. 
In Sec. III we present the equations for determining 
the maximally robust unraveling (MRU) for an atom 



làser model. We concentrate upon continuous Markovian 
unravelings, which give ensembles of Gaussian states, 
and also consider unconstrained Gaussian ensembles. In 
Sec. IV we present the numerical solutions for these equa- 
tions, concentrating on the asymptotic behaviour in the 
limit of large nonlinearity x an d phase noise v. The con- 
cluding Sec. V is a discussion of our results and their 
relation to atom làser coherence, and some suggestions 
for future work. 



II. MAXIMALLY ROBUST UNRAVELINGS 



A. Comparison with Other Approaches 



The idea of robustness has it origins in studies of deco- 
herence and the classical limit p|-f7|,|Ï0l |T5[ . Decoherence 
is the process by which an open quantum system be- 
comes entangled with its environment, thereby causing 
its state to become mixed. However, not all pure states 
decohere with equal rapidity. In particular, Zurek JÏ0| de- 
fined the "preferred states" of open quantum systems as 
those states that remain relatively pure for a long time. 
This idea can be thought of as a "predictability sieve" 
|yj . That is, the preferred states are those for which 
the future dynamics are predictable, in the sense that 
there is some projective question (is the system in some 
particular state?) that is likely to give the result "yes" . 

Our approach, as introduced in Ref. Q and applied 
to resonance fiuorescence by one of us and Brady [16|, is 
to find the maximally robust unraveling. This approach 
shares some similarities with other approaches. It has, 
however, a suite of four distinctive characteristics which 
we enumerate below. 



1. Ensembles of Pure States 



First, we considered not a single pure state, but an en- 
semble of pure states. This is appropriate for situations 
where the open system comes to a mixed equilibrium 
state. The ensemble of pure states that we consider must 
be a representation of that equilibrium mixed state. That 
is, the system has a certain probability of being in one 
of those pure states, as in Eqs. (1.1) and (1.2). Recently, 
Diósi and Kiefer |ÏJ] have also considered ensembles of 
pure states in a similar context. 

Without considering such an ensemble it is necessary 
to put some ad-hoc restriction on the pure states consid- 
ered so that they have some relevance to the actual state 
the system is in at equilibrium. For example, as noted 
above, Gea-Banacloche Q considered only pure states 
having the same mean photon number as the equilibrium 
state of the làser model under consideration. 
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2. Physical Realizability 

Second, in Ref. Q we placed a restriction on the en- 
sembles of pure states that we consider: they must be 
physically realizable. By this we mean that it should be 
possible, without altering the evolution of the system, to 
know that its state at equilibrium is definitely one of the 
pure states in the ensemble, but which pure state cannot 
be predicted beforehand. Diósi and Kiefer JÏJ] have con- 
sidered a similar condition, although they do not make 
the connection with physical realizability and measure- 
ment. In this paper we also consider ensembles without 
the constraint of physical realizability, as it is of interest 
to see how active that constraint is. 

We have considered in detail the issue of physical re- 
alizability of ensembles of pure states in the preceding 
paper ||. Here we merely remind the reader of some 
key points and terminology. An ensemble for a system 
obeying a Markovian master equation is physically real- 
ized by monitoring the baths to which it is coupled. This 
leads to an unraveling [ jÏ7[ of the master equation into a 
stochastic equation for a pure state. In steady, state, the 
pure state will move ergodically within some (perhaps in- 
finite) ensemble of pure states. This is how an unraveling 
defines an ensemble, with the weighting of each member 
being the proportion of time the system spends with that 
state. 



3. Survival Probability 

Third, in Ref. & we defined robustness in terms of the 
fidelity or survival probability of the pure states rather 
than their purity. That is, we consider how close the 
states remain to their original state under the master 
equation evolution, rather than just how close they re- 
main to a pure state. This means that Hamiltonian evo- 
lution alone can affect the robustness of states (whereas 
it does not affect their purity, except in conjunction with 
the irreversible terms). It might be thought that this is 
an undesirable feature. However, as will be shown, using 
the survival probability gives results that accord with the 
usual concept of coherence in làsers. This contrasts with 
the results that are obtained using purity, which we also 
consider at the end of this paper (Sec. V C) 

4- Survival Time 

The final aspect of our work that differs from most 
previous approaches |^,[Ï2| Q is that we quantify the ro- 
bustness by the survival time. (This time was previously 
called the fidelity time in Ref 0). It is the time taken 
for the survival probability to fall below some predefined 
threshold. This is as opposed to considering the rate of 
decay of the survival probability at the initial time. That 
rate is actually identical to half the initial rate of decay of 



the purity, and hence is independent of any Hamiltonian 
terms. It is only by considering the robustness over some 
finitc time that the Hamiltonian terms will contribute. 



B. Unraveling the Master Equation 

In this section, we briefly reiterate the discussion in 
Ref. H on how the master equation is unraveled to yield 
a pure state ensemble. The most general form of the 
Markovian master equation is |Ï8| 

K 

p=-i[H,p]+Y J V[c k ]p = C P , (2.1) 
fc=i 

where for arbitrary operators A and B, 

V[A]B = ABA ] - {A f A,B}/2. (2.2) 

We assume this to have a unique stationary state p ss . It 
can be represented in terms of pure states as 

/O ss = ^p„P„, (2.3) 

n 

where the P n are projection operators and the p n are 
positive weights summing to unity. The (possibly infi- 
nite) set of ordered pairs, 

E = {(P n ,p n ):n=l,2,...}, (2.4) 

we will call an ensemble E of pure states. There are con- 
tinuously infinitcly many ensembles E that represent p ss . 
Our aim is to find the 'best' or 'most natural' represen- 
tation for p ss . 

Our first requirement is that the ensemble be physi- 
cally realizable. This is possible if the environment of 
the system is monitored, lcading to a stochastic quan- 
tum trajectory for the system state. Assuming that the 
initial state of the system is pure, the quantum trajec- 
tory for its projector will be described by the stochastic 
master equation (SME) 

dP = dt[C+U{t)]P. (2.5) 

Here the superoperator U, which we will call an unravel- 
ing, does not affect the average evolution of the system, 
but preserves the idempoteney of P. In the long-time 
limit the system will be in som e pu re state P n , with some 
probability p n such that Eq. ( [2.3] ) is satisfied. Since the 
states and weights will depend on the unraveling U, we 
denote the resultant stationary ensemble by 

£ w = {(PM):n = l,2,...}. (2.6) 

For practical reasons explained in we restrict our 
investigation of the (atom) làser to continuous Markovian 
unravelings (CMUs). As was shown in ||, under a lin- 
earization of the dynamics these lead to Gaussian pure 
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states as the members of the ensembles E u . As men- 
tioned above, we will also consider ensembles, in partic- 
ular Gaussian ensembles, which are not constrained by 
the requirement of physical realizability. This is in order 
to see the importance of this requirement in constraining 
the most natural ensembles. 



2. Comparison with Purity 



As noted in Sec. II A 3 above, it is more common in 



discussions of robustness to use purity rather than sur- 
vival probability. The purity of a state at time t can be 
quantified as 



C. Quantifying the Robustness 

1. Survival Probability 

Imagine that the system has been evolving under a 
particular unraveling U from an initial state at time — oo 
to the stationary ensemble at the present time 0. It will 
then be in the state with probability p". If we now 
cease to monitor the system then the state will no longer 
remain pure, but rather will relax toward p ss under the 
evolution of Eq. (2.1). 

This relaxation to equilibrium will occur at different 
rates for different states. For example, some unravelings 
will tend to collapse the system at t = into a pure state 
that is very fragile, in that the system will not remain in 
that state for very long. In this case the ensemble would 
rapidly become a poor representation of the observer's 
current knowledge about the system. Hence we can say 
that such an ensemble is a 'bad' or 'unnatural' repre- 
sentation of p. Conversely, an unraveling that produces 
robust states would remain an accurate description for a 
relatively long time. We expect such a 'good' or 'natu- 
ral' ensemble to give more intuition about the dynamics 
of the system. The most robust ensemble we intèrpret as 
the 'best' or 'most natural' such ensemble. 

In most of this paper we quantify the robustness of 
a particular state by its survival probability S^(t). 
This is the probability that the system would be found 
(by a hypothetical projective measurement) to be still in 
the state P% at time t. It is given by |Ï9| 



^(í)=Tr[(e«i^) 



(2.10) 



S£(í) = Tr[P^ e "P£] 



(2.7) 



Since we are considering an ensemble E u we must de- 
fine the average survival probability 



S u (t) 



Pn S n(t)- 



(2.8) 



In the limit t — > oo the ensemble-averaged survival prob- 
ability will tend towards the stationary value 



(2.9) 



This is independent of the unraveling U and is a measure 
of the mixedness of p ss . 



The ensemble average of this quantity is also initially 
unity, and approaches Tr[pg S ] as i — > oo. Alternatively, 
the purity could be quantified as the maximum overlap 
of any pure state P n (t) with the evolved mixed state: 



vïit) 



Tr[P n (í) (e ct P%)} 



(2.11) 



For Gaussian states (see Sec. III) these quantities are 
simply related by p%(t) = 2/[l + l/p"(í)]- 

The survival probability has a number of advantages 
over purity. First, we motivated our robustness criterion 
from the desire for E u to remain a good description of 
the system once the unraveling ceases. That is, we wish 
to be able to usefully regard the members of the ensemble 
E u as the states the system is "really" in at steady state. 
This is better quantified by the survival probability be- 
cause the purity effectively takes into account only how 
close the state e £t P „ rem ains to some pure state P n (t) 



[introduced in Eq. ( 2.11 )], not how close it remains to 
the original state P n . An ensemble constructed by con- 
sidering the purity would thus, in general, only remain a 
good description of the system by including the determin- 
istic (but not necessarily unitary) evolution of its mem- 
bers from P n to P n (t) after the unraveling ceases. This 
time evolution would negate the idea that the ensemble 
of states P n is the best representation of the system at 
steady state. 

Another reason for preferring the survival probability 
comes from imagining that the unraveling U continues 
after t = 0. In that case one can calculate a conditional 
survival probability, being the overlap of the pure con- 
ditional state with the pure initial state. The ensemble 
average of this conditional survival probability is sim- 
ply the survival probability S%(t) defined in Eq. (2.7) 
above. Thus the concept of survival probability still ap- 
plies even for the conditional evolution. By contrast, the 
conditional purity of the unraveled state would always be 
unity, and consequently has no re lation to the uncondi- 
tional purity defined in Eq. ( 2.10 ). The latter thus has 
no simple interpretation for the unraveled evolution. 

The final reason for pre ferring survival probability, al- 
ready noted in Sec. II A 3, is that it yields results for the 
atom làser that have a clear and simple physical inter- 
pretation in terms of the coherence of the làser output. 
We will show that this is so in the Discussion section. 

One limit in which quite different results are to be ex- 
pected from using purity rather than survival probability 
is that in which the Hamiltonian part of the dynamics 
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dominates. As will be shown, this limit is highly relevant 
for the atom làser. 

Formally, we split the Liouvillian superoperator C as 



where x is a large parameter and 



A' 



ArrP = y^P[cfc]P- 

k=l 

CrcvP = ~Í[H, p]. 

The reversibility of £ rov implies that 

TrL4£ rev P] = -Tr[B£ ICV A] 



(2.12) 

(2.13) 
(2.14) 

(2.15) 



and so Tr[A£ rcv A] = 0, for arbitrary operators A and B. 

To first order in time, both the survival probability 
and the purity depend only upon the irreversible term: 



s(t) = i + m[pc itl p], 

p(t) =l + 2íTr[P£ irr P]. 



(2.16) 
(2.17) 



For longer times, both expressions will (in general) be 
dominated in the large-x regime by the reversible term, 
but in different ways: 

S(t) ~ 1 + X 2 (t 2 /2)Tv[PCl v Pl (2.18) 
p(t) ~ 1 + xí 2 Tr[P(£ irl X rcv - £ TCV £ in .)P}. (2.19) 

The Hamiltonian term directly affects the survival prob- 
ability, but it affects the purity only in combination with 
the irreversible term. 



3. Survival Time 

The above analysis shows that the difference between 
purity and survival probability only shows up at finite 
times. Thus the best way to characterize robustness is 
to look not at the initial rate of decay of the survival 
probability, but at the time it takes to fall below some 
threshold value A satisfying 



1 > A > Tr[p s 2 



(2.20) 



The ensemble survival time for a particular unraveling 
would then be defined as 



t u =min{í : S u (t) = A}. 



(2.21) 



Note that this time is the first time for which S u (t) = A. 
The survival probability is not necessarily monotonically 
decreasing and in some simple examples there will be 
many solutions to the equation S u (t) = A [lq|. 

A natural choice of A, suggested in Ref. p |, is the max- 
imum eigenvalue of p ss : 



A= hm (TrIO 1 /" (2.22) 

n — >oo 

= max{\j e M : p m Qj = XjQj, Qj = Q 2 }. (2.23) 

This can be shown to satisfy A > Tr[/> 2 S ] as follows. Let 
the eigenvalues of p ss be ordered such that A = Ai > 
A2 > A3 . . .. Then 



Í=2 

<A 2 +^AA, 

J'=2 

= A 2 + A(l - A) = A. 



(2.24) 
(2.25) 
(2.26) 



Here the strict inequality holds unless all eigenvalues of 
p ss are equal. 

In the absence of any monitoring of the bath, the pro- 
jector Qi would be one's best guess for what pure state 
the system is in at steady state. The chance of this guess 
being correct is simply A, which is obviously independent 
of time t. Using this A, the survival time t u could thus 
be interpreted as the time at which the initial state P% 
ceases (on average) to be any better than Qi as an esti- 
mate of which pure state is occupied. In other words, the 
ensemble E u is obsolete at time t u . 

In this paper we do not use this choice for A, for rea- 
sons to be explained later. This brings a certain degree 
of arbitrariness into the analysis. However, as we show, 
the most important and interesting results we obtain are 
independent of the choice of A. 

Having chosen a particular value for A, the survival 
time t w quantifies the robustness of an unraveling U. 
Let the set of all unravelings be denoted J. Then the 
subset of maximally robust unravelings Jm is 



J M = {K G J : t k >r u yU e J} 



(2.27) 



As noted above, in practice it may be necessary to restrict 
the analysis to continuous Markovian unravelings D, and 
the corresponding subset Dm ■ Even if Jm has many el- 
ements IZi, IZ2, these different unravelings may give 
the same ensemble E n = E ni = E^ 2 = ■ ■ ■. In this case 
E n is the most natural ensemble representation of the 
stationary solution of the given master equation. When 
we consider ensembles that are not constrained by the 
condition of physical realizability, we will denote the most 
robust of these by E R . That is, we reserve the calli- 
graphic 1Z to denote a robust unraveling. 



III. MRUS FOR THE (ATOM) LÀSER 

A. The Master Equation 

The master equation we use for the (atom) làser is the 
same as that in the preceding paper H. In the interac- 
tion picture, and measuring time in units of the output 
decay rate, it is 
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p= (^[a^lat]" 1 
-iC^a) 2 ^]. 



-PH + TVXVa]^ 



(3.1) 



The parameters TV and C represent excess phase noise 
and self-interaction energy respectiv ely. This h as t he sta- 
tionary solution expressed in Eqs. (1.1 ) and (1.2), with 
mean boson number fi. 

To make progress on this equation we linearized it 
around a mean field by making the replacement 



« = VM + {x + iy)/2, 



(3.2) 



with x and y Hermitian. The linearized master equation 
has a Gaussian solution with moments 





f^mn - 


= {{x m y mn ) sym ) 


(3.3) 


iven by 








Mio(í) 






(3.4) 


Moi(í) 


= Moi(0)- 


XMio(0)(l - w), 


(3.5) 


Mao(í) 


= M2o(0)w 2 


+ í-w 2 , 


(3.6) 


Mii(*) 


= /in(0)tu 


-x{l + MM2o(0)-2] 








+ ^ 2 [l-M2o(0)]}, 


(3.7) 


M02 (í) 


- M02(0) + 


(2 + i/)í-2xmi(0)(l-«í) 






+ 2x 2 {* 


+ [/i 20 (0)-2](l- W ) 





+ [l-M2o(0)](l- W 2 )/2}, (3.8) 
e _t , x = 4/j.C and f = AjjbN. The long-time 



where 

limit of this is a Wigner function 

W ss (x,y) cx exp(-a; 2 /2) 



(3.9) 



with amplitude quadrature (x) variance of unity and 
phase quadrature (y) variance of infinity This is what is 
expecte d a s the linearized version of the stationary state 



of Eq. (1.1 



The conditions for the output of the làser to be co- 
herent, in the sense of having an atom flux much greater 
than the lincwidth (as conventionally defined) are simply 
stated in terms of the dimensionless self-energy x &nd ex- 
cess phase diffusions v jj) 

3/2 



X < M 
v < fi 2 



(3.10) 
(3.11) 



numbers, (J,u,H20, which play the role of U in Eq. (p^q). 
We do not need H02 because the purity of the unraveled 
states implies that 



2 i 

M20M02 — M11 — L - 



(3.12) 



However, it should be noted that the mapping from U to 
M111M20 is, in general, many-to-one. 

The ensemble can thus be represented as 

E u = {(Pg, íi ,ply):x,y}, (3.13) 

where the second order moments of the pure state P^ç 
are determined by the unraveling U. 

The weighting function is flat for y and for x, is given 

by@ 

p u (x) = [27r(l- M20 )]- 1 / 2 exp{-x 2 /[2(l- M20 )]}. 

(3.14) 

It is convenient to use a new notation for the second 
order moments, 



(3.15) 
(3.16) 
(3.17) 



M02 


M20M02 — 


M11 


Mil 




M20M02 — 


M11 


M20 




M20M02 — 


M11 



For pure states satisfying Eq. (3.12), we have (as in the 
preceding paper 0) 



a 



M02 ; & = M11 ; 7 = M20, 



(3.18) 



The different ensembles are now indexed by the pair /3, 7. 
Not all pairs /3, 7 correspond to physically realizable en- 
sembles. The method for determining which do corre- 
spond to PR ensembles is described in the preceding pa- 
per and the constraints that apply are simply 7 > 
and 

(-2x/3 + 2 + i/)(-27 + 2) - 09 + X7) 2 > 0. (3.19) 



C. Survival Probability 



B. The Unraveled Master Equation 

Under a continuous Markovian unraveling the long- 
time solutions for the linearized stochastic dynamics are 
still Gaussian ||. In fact, the evolution of the second 
order moments M20)Mo2,Mn is deterministic. This mcans 
that for a given unraveling U the stationary ensemble 
will consist of Gaussian pure states all having the same 
second order moments. They are distinguished only by 
their first order moments x, = /zio, y = Moii which there- 
fore take the role of the index n in Eq. (2.6). The differ- 
ent ensembles themselves are indexed by another pair of 



We are interested in the survival probability of the 
states P^y- It is convenient to consider the correspond- 
ing Wigner functions, Wg ^(x, y). Obviously the survival 
probability is independent of y so we will drop this sub- 
script, and set y = for ease of calculation. For Gaussian 
states the Wigner function is a bivariate Gaussian distri- 
bution with the moments fi mn defined above. The state 
with initial moments Mmn (0 ) w ill evolve into a state with 



moments jj. mn {t) given by (3.4-3.8). We will denote the 



Wigner function for the former state Wx(x, y, 0) and that 
for the latter W s (x, y, t). The survival probability of the 
state P%f is given by EG] 
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S s (t) = Ti[P s e ct P s ] = 4tt y dxdyW s (x,y,0)W s (x,y,t) 



4n J dxdyj\f(0)exp 
x Aí(t) exp 



^20(0)^02(0) 



(x~x) 2 (i n (0)(x - x)y y 2 



M2o(0)mo 2 (0)-aíii(0) 2 V 2^(0) M2o(0)m 02 (0) 2 Mo2 (0) 
M2o(í)Mo2(í) / (1-iw) 2 , /in(í)(a;-Sïü)(y + xx(l-iü)) (y + xx(l - w)) 2 



(3.20) 
(3.21) 



AÍ20(í)M02(í) ~ Ptll(í) 2 V 2^Í20(í) 



M20(í)M02(í) 



2/i 02 (í) 



where 



TV = 27TW ^20^02 - M11 



This survival probability should be averaged over all x, weighted by the distribution ( 3.14 ) to get 

S u {t) = J dxS s (t)p u (x). 
Thus S(t) is given by a triple Gaussian integral that evaluates to the following: 



S u {t) = 2 x 



1 (agt ~ /??)/[! + (1 ~ lo)Rt] 
(a + a t ){jo + Jt) - {Po + Pt) 2 



(3.22) 



(3.23) 



(3.24) 



where 



2 „i 



R t = a + a t w + IfitXwz + j t x z 



(q + a t w + j3txz) 2 



[(/?o + /3t)(«o + a tw + /3 t xz) - (a + a t )(Po + f3 t w + jtX z )]' 
(a + a t )[(ao + at)(7o + lt) - (Po + Pt) 2 ] 



(3.25) 



where z = 1 — w and q, / 3, 7 ar e as in Eqs. (3.15)— (3.17), 
and /Lt m „ are as in Eqs. (|3.4|) -(3.8). Note that at t = 
the state is pure, so that «o = 1^02, Po = ^11,70 = M20 
as previously. The survival probability S u (t) is thus a 
function of the initial state parameters 70 and /3o, and 
the dynamical parameters v and X- 



D. The Survival Time 



lim x/Tr[p«] = (2*//) 



-1/2 



(3.29) 



Following the general theory described in Sec. |II C 3j , 
we define the survival time t u as the smallest (in this 
case it will be the only) solution to the equation 



From these expressions it i s evid ent that there would 
be a problem in choosing Eq. ( 3.29| ) for A: it is very close 
to the value for Tr[p 2 s ] = (4-kh)~ 1 / 2 . This means that 
the survival time would be equal to the time by which 
the system has relaxed almost to the equilibrium mixed 
state. In particular, its phase would necessarily be poorly 
defined by this time, which means that the linearization 
of the dynamics that we have been using would not be 
vàlid. 



If instead we start with the solution (3.9) of the lin- 
earized dynamics, we have an even worse situation: 



S u {t u ) = A, 
where A is a constant satisfying 
1 > A > Tr[ Ps 2 s ]. 



(3.26) 



(3.27) 



lim ï/Tr[pl l s } 



0. 



(3.30) 



From the solution (1.1) of the nonlinear dynamics, the 

(3.28) 



lower bound on A is, for (i> 1, 

Tr[p 2 s ] = (47r M )- 1/2 . 
In the same limit, the largest eigenvalue for p ss is 



In this case the survival time would always be infinitc, 
which is not helpful. 

Because of these problems, we have not chosen the 
largest eigenvalue of p ss for A. Instead we have inves- 
tigated the dependence of T n on A for various vàlues, 
namely A = 0.5, 0.2, 0.1, 0.05. As will be shown, the most 
robust ensemble, (that with the largest survival time) is 
substantially independent of A. Unless otherwise stated 
we c hoose A to be the midpoint of the two bounds in 
Eq. (|3.27D, namely 
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A = l/2. (3.31) 



E. Unconstrained Gaussian Ensembles 

Finding the most robust PR ensemble E n consists of 
a searching for the maximum t in the region of space 
allowed by the PR constraint. To determine how impor- 
tant this constraint is in determining E 71 , we also search 
for the maximum r in all of /3-7 space (subject only to 
< 7 < 1). The ensemble picked out by this search we 
will call the most robust unconstrained ensemble and dc- 
note E R . Although we call in unconstrained, it is in fact 
constrained to be of the same form as the ensembles re- 
sulting from a continuous Markovian unravelings. That 
is, it consists of Gaussian states with identical second- 
order moments distinguished only by their mean ampli- 
tude and phase. 



IV. RESULTS 

A. Varying \ with v = 

First we present the results for no excess phase noise 
(y = 0) to see the effect of varying the self-energy pa- 
rameter \. Because our results are numerical, we present 
them mostly in a graphical form. 



1. Evolution at x = and x = 50 

Fig. [ï] shows the evolution of various initially purc 
Gaussia n qu a ntum states under the linearized evolution 
of Eqs. (|3.4| ) — ( 3.8 ) . We represent these states by the one- 
standard-deviation ellipses of the Wigner function. In 
each case we choose the initial mean location of the state 
in phase space to be x = y = 0, and, for the last two 
cases, for y = 0, x — ±-\/3/2 as well. 

The first case in Fig. [Wa) is for v — 0, x — 0, and 
an initial coherent state. The ellipses are plotted for 
t = 0, 3, 10. The middle time is the ensemble-averaged 
survival time for an ensemble of coherent states; that is, 
the time at which the ensemble-averaged survival proba- 
bility S(t) drops to 1/2. For the particular case of the co- 
herent state there is no distinction between the ensemble- 
averaged survival probability and the survival probabil- 
ity of a single coherent state S s (t). That is because the 
x-variance 7 of a coherent state is equal to unity, the 
ensemble-averaged x-variance, so that perforce x = 0. 
Note that the only dynamics in evidence here is phase 
diffusion, causing the y-variance of the state to increase. 
For x = v — 0, the coherent state ensemble is in fact the 
most robust ensemble. This can be verified analytically. 
It is also physically realizable, as shown in the preceding 
paper {§. 



The second case in Fig. |l|(b) is again for an initial 
coherent state but with v = 0, x = 50, plotted for 
t = 0, 0.0678, 0.2. Again the middle time is the survival 
time for the coherent state. Note that it is almost two 
orders of magnitude smaller than the coherent state sur- 
vival time for x = 0. The effect of the large x is to rapidly 
shear the state. This is because the a' a 2 nonlinearity 
amounts to an intensity-dependent frequeney shift. The 
coherent state ensemble E^ a \ however, is not the most 
robust ensemble for x = 50. 

The third case in Fig. ^(c) is the most robust uncon- 
strained ensemble -E R for v — 0,x = 50, as determined 
by the numerical method discussed in Sec. III. Threc 
members, x = 0, ±\/3/2, of this ensemble are displayed. 
Note that the t — state is a highly amplitude-squeezed 
state. In fact it is not purely amplitude-squeezed; the 
rr-y-covariance /3 R = fiu is equal to 0.225. In general, 
the angle 9 between the major axis of the ellipse and the 
y-axis is 



1 2/3 
— aretan — 

2 a 



7 



— aretan — 
2 1 



2/?7 



+ f3 2 - 7 2 



(4.1) 



In the limit of small 7 and p this becomes 9 ~ (3j. In this 
case, with 7 R = 0.100, we have 6> R = 1.2°. This angle 
of rotation is almost too small to make out in the figure. 
It is nevertheless interesting that this slight rotation per- 
sists for all x > 0, and that it is actually in the opposite 
direction to the rotation caused by the shearing. That is, 
as the most robust state evolves it passes through a point 
where the squeezing is purely in the amplitude. Because 
the ï-variance 7 R of the states in this ensemble E K is 
less than unity, the different members of E R have differ- 
ent vàlues of x. The threc initial states we show, with 
x = and x — ±-\/3/2, are typical members of the en- 
semble. The states into which these members of the most 
robust ensemble evolve are plotted for t = 0.100 = r 
(the survival time) and t = 0.2 [as in Fig. 0(b)]. Note 
that the survival time is significantly larger that that for 
the coherent state ensemble in Fig. |l](b). 

The final plot, Fig. |ï|(d), shows typical members of 
the most robust PR ensemble E n . That is, the most 
robust ensemble that can be realized by unraveling the 
master equation. It is very similar to the most robust 
unconstrained ensemble _E R , also being highly amplitude 
squeezed with 7^ = 0.092. The three times at which 
its evolution is plotted are t = 0, t = 0.098 = t 71 , and 
t = 0.2 [as in Figs. 0(b) and |(c)]. Note that the sur- 
vival time r TC is marginally smaller than that for the un- 
constrained ensemble, r R . The principal difference from 
Fig. |ï|(c) is that the x-y-covariance has the opposite sign, 
with f3 n = —0.092. This corresponds to a rotation of 
9 n = —0.48°, a rotation that is accentuated as the evo- 
lution progresses. Again, the initial rotation is almost too 
small to see in the figure, but it is a persistent feature for 
large x- 

From Fig. |ï| it is evident that the evolved states from 
the initial state with x = in the robust cases (c) at 
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t = 0.100 and (d) at t = 0.098 are much closer to the ini- 
tial state than the evolved state in the coherent case (b) 
is at time t — 0.0678. This is despite the fact that all of 
these times are the respective survival times at which the 
survival probability drops to 1/2. However, the evolved 
states from the initial states with x — ± -y/3/2 in cases 

(c) and (d) have a lower overlap with their initial states 
than does the evolved coherent state of case (b). This 
clearly illustrates that the survival probability is neces- 
sarily a property of the whole ensemble of states, not of 
a single member. Figure [ï] also shows that the survival 
probability decays for different reasons in different cases. 
In case (a) it decays because the evolved state becomes 
morè mixed, due to phase diffusion. In case (b) it de- 
cays primarily because the evolved state changes shape 
(shearing) while remaining relatively pure. In cases (c) 
and (d) it decays substantially because the mean position 
of the evolved state moves away from that of the initial 
states in phase space. In Fig. |^ we compare the ensemble- 
averaged survival probability S(t) for the four cases in 
Fig. Note that the time scale for case (a) (x = 0) 
differs from that used for cases (b), (c) and (d) (x = 50). 
For short times the survival probability for the coher- 
ent state ensemble (b) is greater than the survival 
probability for the most robust ensembles E n (c) and E n 

(d) . Indeed, the gradient of the survival probability for 
the coherent state ensemble at t — is much less than 
that of the most robust ensembles. This underlines the 
importance of the survival time, rather than the initial 
rate of decay of survival probability, to quantify robust- 
ness. At short times the survival probability generally de- 
cays lin early, due to irreversible processes, as discussed 



in Sec. IIC2. A coherent state minimizcs this form of 



decoherence, resulting in an almost quadratic behaviour 
of 5 |q) (í) for t < x' 1 = 0.02. This can be under stood 
from the asymptotic analytical expression in Eq. ( 2.18| ) 
for the survival probability for a master equation with a 
large reversible term. This expression only applies for the 
survival probability of a single state, but is applicable to 
a coherent state ensemble because all members are effec- 
tively identical. It need not, and indeed does not, apply 
to the more robust ensembles. In comparison with the 
coherent state ensemble, the most robust ensembles are 
affected more by irreversible evolution at short times but 
less by the interplay of reversible and irreversible terms 
at longer times. 



2. Most Robust Unconstraïned Ensemble for varying x- 



Having looked in detail at x = and x = 50 we now 
present an overview for x ranging from 1 to 10, 000. In 
this section we concentrate upon the most robust uncon- 
strained ensemble. In Fig. |^ we plot the second-order mo- 
ments a R , /3 R , 7 R defining the most robust unconstraincd 
ensemble E R , as a function of x- We also plot the sur- 
vival time r R for this ensemble, and, for comparison, the 



survival time r'"^ for an ensemble consisting of coherent 
states. 

For vàlues of x less than about 7.7, the members of 
the most robust unconstraincd ensemble are close to co- 



herent states, with a R 



7 



R 



1 and fi R < 1. As noted 



above, the states are sheared in the opposíte direction to 
the shearing produced by x- At x ~ 7.7 there is a dis- 
continuí ty in all state parameters. Below this value the 
maximum survival time r lies on the boundary 7=1. 
Above this value, what was previously a local maximum 
at some point 7 < 1 becomes a global maximum, hence 
the jump in the parameters. This is shown by the contour 
plots of r versus 7 and fi in Fig. ^|. 

As x becomes large, all of the curves plotted in Fig. || 
tend to straight lines on the log-log plot. It is thus an 
easy matter to read off the following power laws from the 
gradients of these lines: 

(4.2) 
(4.3) 
(4.4) 

These results clearly show that as x increases, the most 
robust sta tes become increasingly amplitude-squeezed. 
From Eq. (4.1) the sealing law for the rotation angle of 
the squeezed state is 



a R ^ 




fi R r 


- x" 1/3 


„ 7 R r 





X 



(4.5) 



These sealings with x can be understood by consid- 
ering the causes of the d ecay i n t he survival probabil- 
ity from the equations (3.4)-(|3~8|). A typical highly 
amplitude-squeezed state member of the most robust en- 
semble has a mean amp litud e-quadrature fluctuation x of 
order unity. From Eq. (3.5), the mean y-quadrature will 
therefore change in a time t <C 1 by an amount of order 
Xt. This will result in the significant decay of the survival 
probability if the change %í is of order the Standard de- 
viation a 1 / 2 of the y-quadrature for that squeezed state; 
in other words, if t = r where 



(4.6) 



This reduction in overlap due to the motion of the mean 
phase of the states is clearly illustrated in Fig. |ï|(c) for 
the initial states with x = ± -y/3/2. The survival prob- 
ability will also be affected by an in crea se in the phase 
quadrature variance /j, 02 . From Eq. (3.8), the dominant 
terms for short times are Ho2(t) — a = —2xfit + x 2 lt 2 - 
Evidently a positive value of the initial x-y covariance fi 
can, at some time í, cancel the increase in the phase vari- 
ance caused by the nonzero initial amplitude variance 7. 
This effcct will maximize the survival probability if the 
cancellation oceurs at a time of order the survival time 
t. This gives the second condition 
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-'fix' 1 - 



(4.7) 



This effect is most easily seen for the x = initial state 
in Fig. |ï|(c), where the phase variance at the survival 
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time is little changed from its initial value whereas the 
phase variance a short time later is significantly changed. 
Lastly, we consider the effect of mot ion and diffusion 
in the x direction. From Eq. (|3.6| ), the amplitude- 
quadrature variance increases at a rate of order unity. 
It will cause a drop in the survival probability once the 
increase is comparable to the i nitia l amplitudc variance 
7, which is at r ~ 7. From Eq. (3.4) the mean amplitude 
x decays to at rate unity but this will only cause a 
significant drop in S{t) when the decrease in amplitude 
is of the order of the amplitude Standard deviation, that 



is for 



7 



-1/2 



which is much longcr. Thus the third 



condition is just 



T~ 7 ~ar 1 . (4.8) 

Once again, the x = initial state in Fig. |ï|(c) shows 
that there is indeed a significant increase in the ampli- 
tude variance at t equal to the survival time. 

The maximum survival time will clearly be when the 
survival times from the effects above which cause decay 
of the survival probability are comparable. The u niqu c 
s olu tions to the three analytical scaling relations ( |4.6| ) 
(4.8) are the scaling laws found numerically and given in 
equations fl4.2|)-(4.4) above. 



Not only does r K scale in the same way as 7 R , it ac- 
tually asymptotes to 7 R for large %■ This is a conse- 
quence of our choice A = 1/2, as will be shown later. 
In any case, the ensemble-averaged survival time clearly 
decreases with Xi so that the nonlinearity causes a loss 
of robustness in the system even under a maximally ro- 
bust unraveling. However, this loss of robustness is much 
worse for othcr ensembles. For example, the coherent 
state ensemble has a survival time that varies as 

rW-X" 1 , (4-9) 

as shown by the dash-dot-dot curve in Fig. ||. Thus for 
large x the description of the làser steady state in terms 
of the highly amplitude-squeezed states of the most ro- 
bust ensemble is much more useful than the conventional 
coherent state descripti on. 

The scaling in Eq. (4.9) can be easily derived from 
Eq. (3.8). Even more simply, it can in fact be d erived 
from the asymptotic analytical formula in Eq. ( 2.18| ) 



for the survival probability for a master equation with 
a large reversible term. With P a coherent state with 
x = and C rcv p = — i[(x/4)x 2 , p] we find for the solution 
S(r) = 1/2, 



= V8 



X 



(4.10) 



Even the coeficient here is a reasonable approximation, 
as Fig. D shows. 



3. Most Robust Physical Realizable Ensemble for varying \ 

Having examined the most robust unconstrained en- 
semble, we now determine the effect of the physical real- 
izability constraint as x varies from 1 to 10, 000. This is 



shown in Fig. [|. It can be seen from this plot that the 
ensemble parameters differ from those in Fig. [3] for all x- 
That is, the PR constraint is active for all X- Ther e is no 
discontinuity in the parameters, because Eq. (|3.19| ) keeps 
the state away from the maximum of r in /3- 7 space. 
This is illustrated clearly in Fig. |[ where the shaded re- 
gions represent the PR states. It is also clear from Fig. |] 
that, for large Xi fi is effectively constrained to be neg- 
ative, which is why we plot ±/3 rather than just /3 in 
Fig. |[ That is, the shearing is in the direction induced 
by the nonlinearity, rather than in the opposing direc- 
tion as adopted by an unconstrained ensemble. The PR 
ensemble is, not surprisingly, more physically reasonable. 

Despite these differences, the scaling laws for a 71 , \f3 n \, 
7 TC and t 11 are the same for the most robust PR ensem- 
ble E n as for the most robust unconstrained ensemble, 
that is 



l'R 



2/3 



JR 



7 



R 



-1/3 
-2/3 



(4.11) 
(4.12) 
(4.13) 



The scalings for a K , 7 K and r n can be derived using 
the same reasoning as in the preceding case. The scal- 
ing for f3 n arises as follows. For robustness the system 
would likc to have (3 positive, as argued above. The con- 
straint forces it to be negative, which is why E^ is always 
constrained, and is situated on the boundary of the PR 
rcgion in f3~j space. For x large and 7 small, the bound- 
ary of the PR region can be found from Eq. ( 3.19 ) to 
be 



(4.14) 



-(3 = x7 2 /4, 
which here scales as x -1 ^ 3 - 

B. Varying v 



We turn now to the effect of excess phase noise v. Fig- 
ure ^ is an overview of the most robust PR ensemble for 
X = and for v ranging from 1 to 10, 000. The behaviour 
is very simple. For v < 2.3 the most robust states are 
coherent states. As v increases they become increasingly 
squeezed states. For all vàlues of v we have j3 = (which 
is therefore not plotted) , indicating that the most robust 
states are purely amplitude-squeezed. The scaling laws 
derived from this plot are 



■R 



,1/2 



,'R 



7 



■R 



-1/2 
-1/2 



(4.15) 
(4.16) 
(4.17) 



This ensemble is not constrained by the PR constraint 
19). These scaling can again be deduced by argu- 



(i 



ments similar to those in Sec. IVA 2. Unlikc the nonlin- 



ear x term, phase diffusion does not cause motion of the 
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mea n po sition of a typical squeezed state. Rather, from 
Eq. ( |3.8| ) , it simply causes the phase-quadrature variance 
to increase linear ly as vt. The survival probability will 
drop significantly in this time if vt is comparable to the 
original phase variance, a. From the increase in the am- 
plitude variance we get r 



7 



as in Sec. [VA 2. 



The maximum survival time occurs when these two times 
are comparable, giving T n ~ v~ x / 2 and a 71 ~ v 1 / 2 , as 
found numerically. 

The survival time decreases with increasing v, and, 
once again, it asymptotes to j n for large v. For com- 
parison we also plot the survival time t'°^ for a coherent 
state ensemble. This scales as 



<■>) 



(4.18) 



so that for large v the most robust ensemble is much more 
robust than the coherent state ensemble. This scaling 
can be derived from the short time asymptotic analytic 
expression in Eq. (2.16). Since the excess phase diffu- 
sion dominates the evolution for v large we have approx- 
imately 



S(t) ~ 1 + vtTr{PV[x/2]P}. 



(4.19) 



Again, this expression only applies for a single state or 
an ensemble such as the coherent state ensemble where 
all members are effectively identical. In the latter case it 
evaluates simply to 1 — vt/A. 



C. Varying A 

The final parameter we wish to consider varying is 
A, which defines the survival time r by the equation 
S(t) = A. All of the rcsults presented so far were for 
A = 0.5. In Fig. H we show the parameters a n and 
r K for the most robust ensemble as a function of x f° r 
v = and for four vàlues of A. For large \ t ne slope of 
the curves are ind epend ent of A. Thus the scaling laws 
established in Sec. [VA are independent of A. As A de- 



creases, the survival time t 71 increases, because it takes 
longer for the survival probability to decay to that level. 

Decreasing A also causes the phase variance a n to in- 
crease, indicating that the most robust states are more 
highly squeezed. This is not unexpected, since the differ- 
ence between the coherent state ensemble and the most 
robust ensemble is expect cd to b e greater at longer times 
by the argument in Sec. IVAl. However, the relative 



increase in a n is far less than the relative increase in r n . 
In other words, the most robust ensemble is only weakly 
dependent on A. Interestingly, because j 71 ~ l/a n , r y n 
decreases as A decreases, while t 71 increases. Thus the 
asymptotic result 7 TC ~ r TC can only be true at one value 
of A, namely A = 1/2. 

Figure ^ presents the same information as Fig. |^ but 
for x = and varying v and A. Once again the scaling 
laws established in Sec. IV B are found to be independent 



of A, and in this case the different vàlues for a n appear 
to asymptote. In this case, the value for v above which 
the coherent state ensemble ceases to be the most robust 
ensemble increases for decreasing A. Above these vàlues 
of v the amplitude-squeezing in the most robust ensem- 
ble is always decreased as A is decreased. However, the 
difference is small (and may vanish as v — ► oo), so that 
the equation t 71 ~ j 71 is again vàlid only for A = 1/2. 
The sum of these results justifies our use of the single 
value A — 1/2 for most of this work. 



V. DISCUSSION 
A. Summary 

The atom làser is an open quantum system with rich 
dynamics. In this paper we have explored a new way 
of characterizing those dynamics: finding the maximally 
robust unraveling || . This yields the most robust physi- 
cally realizable ensemble E n of pure states P n that sur- 
vive the best. By "surviving", we mean remaining un- 
affected by the system dynamics. This ensemble is, we 
have argued, the most natural representation of the sta- 
tionary state matrix p ss of the làser; if one wished to 
regard the làser as being "really" in a pure state, then 
the most natural states to choose are the members of this 
ensemble. Although it is a time-independent ensemble, 
it is drastically affected by alterations in the dynamics of 
the atom làser that do not change the stationary state 
matrix. 

We considered a simple model for the atom làser in 
which p ss is a Poissonian mixture of number states of 
mean /i. Working in the linearized regime, we identified 
two relevant dynamical parameters that may be varied 
without altering this stationary state. The first is x, 
which is proportional to the strength of sclf-interaction 
of the àtoms in the làser. The second is v, which is pro- 
portional to the excess phase diffusion of the làser above 
the Standard quantum limit. 

For x — arL d v small, the most robust ensemble 
was found to consist of coherent states, with mean boson 
number \x but with all possible phases. This is the most 
common representation of the state of an optical làser, 
and so it is not surprising. In terms of the parameters we 
used in the paper, the ensemble consists of Gaussian pure 
states with phase quadrature variance a = 1, amplitude- 
quadrature variance 7 = 1, and amplitude-phase covari- 
ance ,9 = 0. 

As the self-energy x is increased the most robust states 
cease to be coherent states. In fact, for any nonzero 
value of x, not only are the coherent states not the most 
robust state; in addition they are not even physically 
realizable ||. For large vàlues of x thfi most robust 
states P n are very highly amplitude-squeezed states with 
amplitude-quadrature variance 7^ scaling as x~ 2 ^ 3 an d 
phase quadrature variance a n scaling as x 2 ^ 3 . The same 
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effect occurs for large vàlues of v, with scalings of //~ 1//2 
and v 1 !" 1 respectively. 

It is not known what value of v would be appropriate 
to model a realistic àtom làser. However, it was argued 
in Ref. H that a typical value for \ might be 1000. This 
implies that the most natural description of an atom làser 
would be in terms of higl·ily amplitude squeezed stated, 
with the Standard deviation in the amplitude quadrature 
being of order 0.1. Excess phase noise would only in- 
crease the amount of squeezing in the states in the most 
robust ensemble. 

As noted above, our analysis was based on a linearized 
approximation for the làser dynamics. This is only vàlid 
if the states under consideration have a well-defined co- 
herent amplitude. As x or v are increased indefinitely 
and the most robust states become more amplitude- 
squeezed, this approximation will clearly break down. 
Specifically, it will break down when the phase variance 
predicted by the linearized analysis is of order unity; that 
is, when the the phase quadrature variance is of order 
the mean boson number fi. From the above scalings, for 
the linearization to remain vàlid we requirc 



X < M 
v < /i 2 



3/2 



(5.1) 
(5.2) 



Although wc cannot say with confidence what the most 
robust states are when the linearization breaks down, we 
do know that they must be states without a well-defined 
coherent amplitude (because that is why the linearization 
brea ks down). Therefore the conditions in Eqs. (|0]) and 
( |5.2| ) also represent the conditions for the most robust 
states to be states with well-defined coherent amplitudes. 
In other words, if and only if these conditions are satis- 
fied, the most natural description of the atom làser is in 
terms of states with a mean ficld. 



B. Interpretation 

We can now finally state t he m ost important result of 
this paper. The conditions ( |5.l| ) an d (|5.2] ) arc iden tieal 
to the previously stated conditions Q3.10] ) and (3.11) for 

Here we mean 



the output of the device to be coherent, 
coherent in the sense that the output is quantum de- 
generate, with many bosons being emitted per coherence 
time. Without this condition the device could not be 
considered a làser at all, as its output would consist of 
independent àtoms rather than a matter wave. 

The significance of this result is that there is a perfect 
correspondence between the 'best ' pure states for describ- 
ing the làser, and the coherence of its output. If the most 
robust states have a well-defined coherent amplitude, like 
coherent states, then the output is coherent. If the most 
robust states do not have a well-defined coherent ampli- 
tude, like number states, then the output is not coherent. 
This profound result establishes the usefulness of maxi- 



mally robust unravelings as an investigational tool for 
open quantum systems. 

It must be emphasized that the link between the pres- 
ence or absence of a mean field inside the làser, and the 
presence or absence of quantum coherence in the làser 
output, is not due to any simple relationship of defi- 
nitions. Finding the maximally robust ensemble is, as 
the diligent reader will appreciate, a very involved pro- 
cess completely different from calculating the first-order 
coherence function. In particular, the average survival 
time for the members of the most robust ensemble has in 
general no relationship with the coherence time. 



C. Comparison with Purity 

It is worth pointing out that the relationship we have 
established between robust mean-field states and quan- 
tum degeneracy would not have been found had we used 
purity rather than survival probability as the basis of our 
definition for the most robust ensemble. Although there 
are no great differences between the two definitions as one 
varies v, there is a great difference as one va ries x- This 

as x 



is to be expected from the analysis in Sec. [I C 2 



seales the self-energy Hamiltonian, whereas v represents 
irreversible phase diffusion. 

To prové this point we have calculated the ensemble 
that maximizes the time it takes for the avera ge purity 
of the member states [as defined in Eq. (pTlOb l to drop 
to 1/2 under the master equation evolution. We plot 
the parameters for this ensemble as a function of x i n 
Fig. |^. For comparison we also plot the phase quadra- 
ture variance a TC and the survival time of the most 
robust ensemble as previously defined, in terms of sur- 
vival probability. The ensemble parameters when we use 
purity obey sealing laws for large x, but they are different 
from those sealing laws obtained when using the survival 
probability (Sec. V A 2): 



JR. 



7 



TZ' 



x 1/2 , 

-1/4, 

x- 1/2 . 



(5.3) 
(5.4) 
(5.5) 



As expected from Sec. II C 2, the purity half-life is much 



longer than the survival time for large x- Here we use TZ' 
rather than TZ to emphasize that we are using a different 
mcasurc of robustness. 



The scalings in Eqs. (5.3)— (5.5) can be derived ana- 
lytically. For Gaussian states with moments /Lt m n(í), the 
purity at time t is given by 

Tr [p 2 {t)\ =p(t) = [M2o(í)Mo 2 (í) - A*nWr 1/2 - (5-6) 
For v = 0, 7 •< 1, [3 ~ 1, y_3> 1 and t -C 1, as appro- 



priate here, the solutions ( |3 . 6|) — (|3 . 8| ) , together with the 
condition p(t) = 1/2, yield the following equation for r 

3 « 2(1 + /3 2 )r/ 7 - 2 X fir 2 + 2 X 2 7 r 3 /3 + *V/3. (5.7) 
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It is clear from the term 0(r 4 ) that r will scale as x 1 ^ 2 · 
To maximize r, the terms 0(t) and 0(r 3 ) imp ly that 
7 should scale as x" 1 ^ 2 in accord with Eq. ( |5.5| ) . The 
terms O(r) and 0(t 2 ) then imply that (3 should be pos- 
itive, and of order unity. Indeed, for the unconstrained 
Gaussian ensemble E R we find (3 » 1.8. With the con- 



straint of Eq. ( 4.14 ), we get (3 negative and of order unity, 
as stated above. 

The condition for the best purity-preserving states to 
have a well-defined coherent amplitude is oP^ <C /i, which 
from Eq. (|5.3|) gives 



(5.8) 



This implics that there is a range of interaction strengths 
M 3 ^ 2 ~ X *C fJ 2 for which the purity analysis delivers a de- 
scription of the làser in terms of states with a mean field 
even though the làser output is no longer coherent in the 
sense defined above. This regime can be interpreted in 
terms of a the non-standard concept of conditional coher- 
ence^ explored in detail in the preceding paper ||. The 
bàsic idea is well illustrated by Fig. If one knows the 
mean amplitude of the state with an uncertainty much 
less than unity, as in Fig. |ï|(d), then the direction that 
it will move in phase space can be predicted with accu- 
racy. This motion (which amounts to different freqüèn- 
cies) can then be taken into account in experiments the 
output. Thus the spread in freqüències due to spread in 
amplitude can be compensated for (up to a point). 



D. Comparison with Quantum State Diffusion 

A particular PR ensemble of interest is that gener- 
ated by the unraveling known as quantum state diffusion 
(QSD) |pï| , p2| . This is merely a particularly simple and 
natural type of continuous Markovian unraveling. It has 
been suggested |Q that the corresponding ensemble is 
a good candidate for the most robust ensemble. We in- 
vestigated this ensemble in the preceding paper || and 
found analytically that its parameters (3 and j have ex- 
actly the same scaling as the PR ensemble E n based on 
maximizing the robustness as measured by purity. That 
is, with x, 



and with ía 



a QSD 
^QSD 
7 QSD 



vQSD 



-1, 

V2X- 1/2 - 

7T 1/2 ' 
o, 



^qsr 
7 QSD „ 72^~ 1/2 . 



(5.9) 
(5.10) 
(5.11) 



(5.12) 

(5.13) 
(5.14) 



Consequently, the QSD ensemble E^ SB scales with x 
quite differently from the maximally robust ensemble E^ 



according to our definition based on maximizing the sur- 
vival time. Thus unlike E n , but like E n , the coherence 
of its members does not have a direct correspondence 
with the làser output coherence (in the conventional, un- 
conditional, sense). 

The correspondence (at least in scaling laws) between 
j^QSD anc j j^R - 1S ac tually in contrast to the result found 
by Diósi and Kiefer (for a different system) They 
found that PR states minimizing the loss of purity were 
different from states produced by QSD. However, as 
noted earlier, they considered only the initial rate of loss 
of purity, which is insensitive to Hamiltonian terms. If 
they had considered maximizing the half-life of the pu- 
rity, as we have, they may have obtained a different re- 
sult. 



E. Future Work 

There are at least three future directions for this work. 
First, the insights into the atom làser that the maximally 
robust unravelings analysis offers suggests that this tech- 
nique could be applied fruitfully to other open quantum 
systems. It has already been applied to fluorescent àtoms 
|16| , and could also be applied to other quantum optical 
systems |ï^] , and other models for Bose-Einstein conden- 
sates in equilibrium with a reservoir (23). These are all 
systems with nontrivial dynamics, which could be more 
fully appreciated by determining the maximally robust 
unraveling. 

Second, the difference between the analyses based on 
survival probability and purity deserves further investi- 
gation. As we showed, the purity analysis gives a de- 
scription of the làser mode in terms of states with a well- 
defined coherent amplitude for high vàlues of x where the 
survival analysis does not, and where the output is not 
coherent in the conventional sense. Nevertheless the re- 
sults do make sense in terms of conditional coherence . 
Perhaps it is because purity is unaffected by the motion 
of the mean position of the states in phase space that it 
refleets conditional coherence, which relics on knowledge 
of that motion to define the output mode. By contrast, 
the survival probability is affected by the motion of the 
states, and hence refleets conventional coherence that av- 
erages over the different freqüències of rotation. 

Finally, there are other approaches to quantifying the 
robustness of unravelings apart from the survival proba- 
bility and the purity. For example, one could measure 
how quickly the unraveling purifies the state, or how 
sensitive the purity is to imperfections in the unravel- 
ings. Related ideas have recently been explored ||Ï5| . p4| . 
These ideas could be best investigated in systems some- 
what simpler than the atom làser we have considered 
here. This would give an indication for the robustness of 
the idea of robustness; that is, how sensitive the maxi- 
mally robust unraveling is to the definition of robustness 
used, and which definitions agree. 
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To conclude, the clear and simple interpretation for [24] Z. Brady, B.Sc. Honours thesis (Grifhth University, 
the results we have obtained here for the atom làser viri- 2000). 
dicates our conviction Q that maximally robust unravel- 
ings will have an increasing rolc as a tool for understand- 
ing the dynamics of open quantum systems. 
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FIG. 1. The evolution of (initially pure) Gaussian quan- 
tum states under the linearized làser master equation for four 
different cases. The states are represented by the one stan- 
dard-deviation ellipse of the Wigner function. In the all cases 
we choose the initial mean location of the state in phase space 
to be x = y — 0, and for the last two we additionally have 
x — ±-\/3/2. For all four cases the excess phase diffusion is 
v — 0. For case (a) we have \ — and an initially coherent 
state (which forms the most robust ensemble in this case). 
For case (b) we have x = 50 and again an initially coherent 
state. For case (c) we have x = 50 but the initial states are 
members of the most robust unconstrained ensemble E R for 
this X- F° r case (d) we have x = 50 but the initial states are 
members of the most robust PR-constrained ensemble E 71 for 
this x- In a U cases the black ellipses are for t = 0, the dark 
grey ellipses for t — r (the appropriate ensemble-averaged 
survival time), and the light grey ellipses for a still later time. 
Details of these times are given in the main text. 
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FIG. 2. The decay of the ensemble-averaged survival prob- 
ability in time for the four cases represented in Fig. |ï|. The 
horizontal axis measures time í. For case (a) it is scaled in 
units of the bare lifetime of the làser mode, and for cases (b), 
(c), and (d) it is scaled in units 100 times smaller. That is, 
the survival probabilities actually drop much more quickly for 
the last three cases. 
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FIG. 3. The parameters for the most robust unconstrained 
Gaussian ensemble E R as a function of x with v — 0. These 
parameters are the phase quadrature variance q r (dotted 
line), the amplitude-quadrature variance 7 R (dashed line), 
the covariance (3 R (dash-dot line) and the survival time r R 
for the members of this ensemble. For comparison, we also 
plot the survival time r' a ^ (dash-dot-dot line) of a coherent 
state ensemble. Both survival times are in units of the bare 
lifetime of the làser mode. 
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FIG. 4. Contour plots of the survival time r as a function 
of 7 and f3. In (a) v — and \ ~ 7-7 and in (b) v = and 
X — 50. In each plot the heavy curves represent contours of 
t (in units of the bare lifetime of the làser mode) and the 
shaded region represents states that are physically realizable 
(PR). Crosses mark the positions of the maximally robust 
ensembles (MRE). 
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FIG. 5. The parameters for the ensemble E arising from 
the maximally robust unraveling 1Z as a function of \ with 
v = 0. As in Fig. ^ we plot a 71 (dotted line), 7^ (dashed line) 
and ±[3 K (dash-dot lines). We also plot the survival time t 71 
(sòlid line) of this ensemble and, for comparison, the survival 
time t'°^ (dash-dot-dot line) of a coherent state ensemble. 
Both of these times are in units of the bare lifetime of the 
làser mode. 
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FIG. 6. The parameters for the ensemble E n arising from 
the maximally robust unraveling K as a function of v with 
X = 0. As in Fig. U we plot a n (dotted line), "/ n (dashed 
line), and the survival time r 71 (sòlid line). We do not plot 
(3 n because it is identically zero. For comparison we also plot 
the survival time r' Q ^ (dash-dot-dotted) of a coherent state 
ensemble. Both of these times are in units of the bare lifetime 
of the làser mode. 
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FIG. 7. The parameters for the ensemble E arising from 
the maximally robust unraveling 1Z as a function of \ with 
v = and for various A. The rising lines are a 71 and the 
falling lines are r 71 (in units of the bare lifetime of the làser 
mode). The vàlues of A are 0.5 (sòlid line), 0.2 (dashed line), 
0.1 (dash-dot line), and 0.05 (dotted line). 
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FIG. 8. The parameters for the ensemble E K arising from 
the maximally robust unraveling K as a function of v with 
X = and for various A. The rising lines are a n and the 
falling lines are T n (in units of the bare lifetime of the làser 
mode). The vàlues of A are 0.5 (sòlid line), 0.2 (dashed line), 
0.1 (dash-dot line), and 0.05 (dotted line). 
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FIG. 9. Parameters for the maximally robust ensemble for 
v — as a function of x as i n Fig. ^ but using purity as a 
measure of robustness. As in previous figures we plot a 71 
(dotted line), r y TÍ (dashed line), /3 TC (dash-dot line) and t 71 
(sòlid line). Also shown for comparison are the ot 1 (rising) 
and r n (falling) curves from Fig. |5| as dash-dot-dotted curves. 
Both times are in units of the bare lifetime of the làser mode. 
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